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"To divide is to leave something undivided. 
To discriminate between alternatives is to 
leave something which is neither alterna- 
tive." 

Chuang-tzu 
Abstract: 



These lecture notes form a primer on the theory of cosmic microwave back- 
ground (CMB) anisotropy formation. With emphasis on conceptual as- 
pects rather than technical issues, we examine the physical foundations of 
anisotropy evolution in relativistic kinetic and perturbation theory as well as 
the manifestation of these principles in primary and secondary anisotropies. 
We discuss gauge choice and gauge invariance and their use in understand- 
ing the CMB. Acoustic, gravitational redshift and ionization effects have 
robust signatures in the CMB spectrum and may allow determination of 
classical cosmological parameters as well as reveal general distinctions be- 
tween models for structure formation. We develop the tight and weak cou- 
pling approximations as analytic tools to help understand these effects and 
the robustness of their signatures. 
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1 Introduction 

More than three years has passed since cosmic microwave background 
(CMB) anisotropics were first detected at large angular scales in the COBE 
DMR sky maps (Smoot et al. 1992). Since then much progress has been 
made in probing the anisotropy spectrum at smaller scales (see e.g. the 
compilations by Bond 1995; Scott, Silk & White 1995; Steinhardt 1995). 
Satellite missions now being proposed would be able to map the anisotropy 
definitively down to a fraction of a degree. To exploit these recent and po- 
tential experimental advances, we need an accurate and accessible theory of 
anisotropy formation. The first part of this task was largely accomplished in 
the 1980's with the ground-breaking work of Wilson & Silk (1981), Vittorio 
& Silk (1984), and Bond & Efstathiou (1984). Detailed numerical calcula- 
tions in perturbation theory allow accurate predictions of the anisotropy in 
most models for structure formation. Summaries of the current state of the 
art in numerical codes can be found in Bond (1995), Hu et al. (1995) and 
Ma & Bertschinger (1995). 

Numerical work however is notoriously inaccessible to the uninitiated. 
Moreover, it obscures the true potential of CMB anisotropy spectrum for 
cosmology by limiting itself to specific models. Indeed, CMB anisotropics 
are sensitive to classical cosmology parameters such as the matter and radi- 
ation content, baryon fraction, expansion rate, curvature, and cosmological 
constant (Bond et al. 1994; Seljak 1994; Hu & Sugiyama 1995a) as well as 
the model for structure formation (Hu & Sugiyama 1995b; Albrecht et al. 
1995, Crittenden & Turok 1995; Durrer et al. 1995; Magueijo et al. 1995). 
In these lecture notes, we will review the basic concepts of anisotropy for- 
mation and their consequence for cosmology. 

We begin in §2 with relativistic kinetic and perturbation theory which 
give the basic physical laws that govern anisotropy formation. Pure reduc- 
tionism of this sort however would miss the true use of CMB anisotropics. 
In §3 and §4, wc map out the manifestations of the basic principles be- 
fore and after recombination, i.e. the primary and secondary anisotropics. 
These anisotropics possess information on both the cosmological model and 
structure formation. From this decomposition and reconstruction, the re- 
sulting complex structure of anisotropics can be understood in a general 
model-independent manner. 
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2 Physical Foundation 

Two ingredients are necessary to describe the evolution of cosmic microwave 
background anisotropies: relativistic kinetic and perturbation theory. Ki- 
netic theory describes the radiation transport properties of the CMB pho- 
tons in the metric perturbed by gravitationaUy unstable density fluctua- 
tions. These fluctuations are evolved through relativistic perturbation the- 
ory. Combined, the two considerations yield a complete system that is re- 
sponsible for anisotropy formation in any model where structure formation 
proceeds by gravitational instability. 

2.1 Relativistic Kinetic Theory 

Conceptually, relativistic kinetic theory is identical to the familiar non- 
relativistic case: the phase space distribution is conserved along geodesies 
save for a collision term due to scattering, 

^ ^ ^ < < ^tll = r\f] (^) 

dt ~ dt dx' dt dpdt dr dt ^ ' 

Here 7^ are the direction cosines of the photon momentum p. Before a red- 
shift of ~ 1000, CMB photons were hot enough to ionize hydrogen. 
Consequently, the dominant interaction process for the CMB at early times 
was Compton scattering off free electrons. Due to the higher electron den- 
sity in the early universe, the Compton mean free path was quite short, 
much smaller than the particle horizon at that time. The horizon at z^. sub- 
tends a degree or less for Qq-\-Qa ^ 1- Thus on scales relevant for observable 
anisotropies, the photons were tightly coupled to the electrons, which in turn 
were tightly coupled to the protons by Coulomb interactions. At 2;*, neutral 
hydrogen formed through "recombination" and the photons last scattered. 
Unless the universe suffered reionization at high redshift, fluctuations in the 
CMB at recombination were frozen in at z^ and await observation today. 
There are several features of Compton scattering worth noting 

1. Scattering couples the photons to the baryons and forces perturbations 
in their number and hence energy density to evolve together. 

2. Scattering isotropizes the photons in the electron rest frame thus cou- 
pling the local CMB dipole to the electron velocity. 

3. In the Thomson limit, there is no energy transfer in scattering. Energy 
exchange only occurs to i.e. 0{Te/me). 

4. There is no change in photon number through Compton scattering. 

Scattering thus governs the intrinsic temperature perturbations at last scat- 
tering and the dipole or bulk velocity perturbation. Since it does not change 
the net energy or photon number in the CMB to lowest order, spectral dis- 
tortions to the blackbody do not arise in linear theory. Unless the electrons 
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have been heated significantly above the temperature of the CMB or pho- 
tons and/or energy has been dumped into the CMB from an external source, 
the spectral information can be ignored. 

The remaining subtlety is that the photons propagate in a space-time 
that is distorted by density fluctuations: we must employ the geodesic equa- 
tion in the presence of perturbations. This leads to gravitational redshift 
effects from the dp/dt term in equation (1) which can also generate fluctu- 
ations in the CMB. The d^/dt term represents gravitational lensing and 
gives a first order contribution only in a curved universe (see §3.4). By inte- 
grating equation (1) over frequencies, one obtains the Boltzmann equation 
for the evolution of temperature perturbations 0(r7,x, 7) = AT/T, 



d d .jd 



dr] ' " dx'^ ' ' J ^ '^'^ ^ ^ 

where So and Sc are the sources from the gravitational redshift and 
Compton scattering and overdots are derivatives with respect to confor- 
mal time r] = J dt/a. For a derivation from first principles of the Comp- 
ton collision term Sc, see Hu, Scott & Silk (1994), Dodelson & Jubas 
(1995), and Kosowsky (1995). If the angular dependence and polarization 
of Compton scattering is ignored it reduces to f[Oo — O — ^ivl\. Here 6*0 is 
the isotropic temperature fluctuation and the differential Compton optical 
depth f = XeUeCrTCt with Xe as the ionization fraction, Ue as the electron 
density and (Jt as the Thomson cross section. Implicit in So and Sc are 
the three fundamental sources of anisotropies in the CMB: 

1. Gravitational redshifts from the presence and evolution of metric fluc- 
tuations. 

2. Hot and cold spots from the intrinsic temperature at last scattering. 

3. The Doppler effect due to the velocity of the last scatterers. 

Of course, this high level description is not very practical. Gravitational 
instability controls the evolution of density, velocity, and metric perturba- 
tions in the universe. All of these effects are thus related by relativistic 
perturbation theory. 



2.2 Relativistic Perturbation Theory 

General relativity tells us that matter moves in a space-time perturbed by 
fluctuations in the matter density itself. Thus the Einstein equations reduce 
conceptually into two pieces. The stress-energy tensor of the total matter is 
covariantly conserved in the perturbed metric, T'^'^.^ = 0, and matter fluc- 
tuations are the source of metric perturbations via a generalized Poisson 
equation. For the former, 1^ = gives number or energy density conserva- 
tion, i.e. the continuity equation. The spatial components give momentum 
conservation, i.e. the Euler equation. 
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Fig. 1. Gauge ambiguity. Due to the change in density with the expansion, even a 
homogeneous FRW universe can appear to have a density perturbation if the time 
shcing is warped. The synchronous and total matter gauges choose the time slicing 
to correspond to the rest frame of the freely falling and total matter respectively. 
The Newtonian gauge chooses slicing such that the expansion rate is shear free. 

2.2.1 Gauge Choice 

To explicitly implement these principles, we need to treat a subtlety due 
to gauge freedom in general relativistic perturbation theory. In order to 
define a perturbation, we must specify the relation between the physical 
spacetime and the hypothetical unperturbed background. The difference 
between quantities at the same coordinate values is deemed a perturbation. 
Since this choice is not unique, one must generally fix a gauge before making 
any calculations or interpretations regarding perturbation evolution. 

Practically speaking, a gauge choice involves fixing the constant-time hy- 
persurfaces and the spatial grid on these surfaces. Suppose we warp the time 
slicing from some fiducial choice by ry = rj + Te^^'^ (see Fig. 1). The density 
of a fiuid of particle X on this slicing becomes px = Px — pxTe'^^'^ for the 
same coordinate value. Because the expansion makes densities decrease in 
time as 

Px = -3(px +Px)-, 
a 

in the new frame the density fluctuation is 

~Sx = Sx+ 3{l+px/px)-Te^-''. (3) 

a 

As a simple example of gauge ambiguity, consider a Friedmann-Robertson- 
Walker universe with no perturbations 6x = 0. Even in this density 
fluctuation appears in the warped time slicing. One might argue that this 
involves a silly choice of time slicing and that observers may set up a sensible 
coordinate system that removes ambiguities by exchanging light signals. 
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This is true. All "sensible" choices of coordinate systems agree on the nature 
and value of the density perturbation within the horizon by requiring that 
the warping satisfy drj / rj = T / rj <^ 5 il krj ^ 1. However outside the horizon, 
there is no "sensible" way to set up the coordinate system. In particular, 
different choices of time slicing can produce seemingly different results for 
perturbations and their evolution in this regime. 

There are two basic approaches to fixing a gauge or coordinate system. 
We can choose a set of preferred observers and set the coordinate frame 
to be their rest frame, i.e. employ "Lagrangian" coordinates. The popular 
and useful synchronous gauge is one such representation. Here we choose 
freely falling observers to define the frame, e.g. the coUisionless cold dark 
matter. This rest frame analysis is convenient in that it simplifies the evo- 
lution equations and their solution. An extension of this idea is employed in 
the total m,a,tter gauge. Here the coordinate system coincides with the rest 
frame of the combined relativistic and non-relativistic matter. This choice is 
computationally useful in the early universe where radiation dominates and 
does not necessarily follow the perturbations in the non-relativistic matter. 

The computational convenience of these "Lagrangian" techniques is 
offset by difficulties in obtaining physical intuition for processes such as 
anisotropy formation. Perturbations grow by dilation effects due to distor- 
tions of the coordinate grid rather than simple causal mechanisms such as 
gravitational infall and redshift. An alternate choice is to take an "Eule- 
rian" coordinate system. More specifically, instead of employing a set of 
preferred observers, fix the coordinate system by a geometric property. One 
useful choice is the Newtonian gauge. The coordinate frame here corre- 
sponds to zero-shear hypersurfaces where the expansion appears isotropic. 
Inside the horizon, this analysis reduces to the standard non-relativistic 
Newtonian treatment. Thus, our physical intuition from the non-relativistic 
theory most easily carries over to a complete treatment in this gauge. 

Since there is a one-to-one mapping between gauges whose coordinates 
are entirely fixed, these treatments are entirely equivalent and give the same 
result for any observable quantity, e.g. CMB anisotropics. Ambiguities only 
arise if the coordinate system is not entirely fixed. In this case, the one- 
to-many mapping produces unphysical gauge modes which historically has 
produced much confusion (see Press & Vishniac 1980). Gauge modes can 
be avoided by completely fixing the coordinates. The definitions of the total 
matter and Newtonian gauges already do so. On the other hand, a syn- 
chronous coordinate system, defined by a set of freely falling observers, 
requires us to specify which set through the initial conditions. In particular, 
if the initial velocity of the observers is fixed then there are no gauge modes. 
In practice, one usually takes the rest frame of the CDM. 



Concepts in CMB Anisotropy Formation 7 

2.2.2 Gauge Invar iance 

Due to the computational convenience of the "Lagrangian" gauges and the 
intuitive nature of the "Eulerian" ones, it is often useful to switch between 
gauges in the midst of a calculation or even employ different gauge rep- 
resentations in a single set of evolution equations. A useful technique for 
this purpose is provided by the gauge invariant method (Bardeen 1980; Ko- 
dama Sz Sasaki 1984; Mukhanov, Feldman & Brandenberger 1992). It is 
called gauge invariant not because a density perturbation, say, remains a 
density perturbation under a gauge transformation. Rather, a density per- 
turbation in a specific gauge is represented in a coordinate free manner, 
i.e. in a manner valid in an arbitrary gauge. In practical terms, this just 
means that whatever gauge is chosen, one has a systematic way to mix in 
quantities that represent perturbations in a different gauge. 

Let us make these considerations more concrete. The most general form 
of a metric perturbed by scalar fluctuations a given Fourier mode k is 
(Bardeen 1980) 

^oo = -a2[l + 2A^e*-"], 
goj = a^B^ikje^^'', 

9ij = a''{5ij + [2H^5ij + 2H^{-kikj + 5i,/3)]e^'^-"}, 

where for illustration purposes, we have taken a flat geometry and super- 
script G is meant to remind the reader that these quantities are gauge 
dependent. By fixing the gauge, we can eliminate two out of the four metric 
terms. In the Newtonian gauge, = = 0, and it is traditional to 
label the Newtonian potential = ^ and the space curvature perturba- 
tion = <P. Synchronous gauge requires = = 0, whereas the total 
matter gauge takes B^ — and H'^ = where is the velocity of 
the total matter. Ordinarily, one now writes down the Einstein equations 
G^jj = SnGT^^y which contain the conservation T^'^.^ = and "Poisson" 
equations, with a specific choice of gauge. For a fiuid of particle X, the 
stress energy momentum tensor is 

n--(l + 5fe^'^-"), 

T° =(px+Px)Vf M.)e^'^•^ (5) 
T'j = {px + 5p^e'''-^)5'j+pxnx{-k% + 5^)e*•'', 

where Ux is the anisotropic stress of the fluid. 

As an example of the gauge invariant program, let us write the total 
matter gauge density perturbation in a gauge-independent way. The total 
matter gauge condition is satisfied by shifting the time slicing from an ar- 
bitrary gauge by T = {V^ - B'^)/k [Kodama & Sasaki 1984, eq. (3.3a); 
Hu 1995, eq. (4.88)]. Equation (3) tells us that the density perturbation 
transforms as 



(4) 
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Sx = Sx + Hl+Px/px)-{V^ - B'')/k. (6) 

In an arbitrary gauge G, this represents what the density perturbation would 
be in the total matter rest frame. This "gauge invariant" definition can be 
used to employ total matter gauge variables in a Newtonian or synchronous 
gauge treatment: 

5l^5^ + 2>-{l+px/px)V^/k 

« (7) 

= 4 + 3-(l+px/px)^T^A. 

a 

If fiuctuations are adiabatic, the number density fiuctuations of the matter 
and radiation evolve together, i.e. their bulk velocities are equal. In this 
case, Vrp is equal to the freely-falling cold dark matter velocity, commonly 
defined to be zero in synchronous gauge. Equation (7) then tells us that 
the total matter gauge and synchronous gauge density perturbation are nu- 
merically equivalent. In a matter dominated (pressureless) universe, there 
is no fundamental scale in these rest frame evolution equations and density 
fluctuations evolve in a scale free manner. This is not true for the Newto- 
nian gauge. Here fluctuations grow by the gravitational infall of the matter 
into potential wells. Due to causality, this introduces a fundamental scale, 
the horizon scale, into the evolution equations. For example in the case of 
adiabatic fluctuations, density perturbations are constant outside the hori- 
zon and only grow (or decay) after the horizon has grown larger than the 
wavelength. 

Now let us see how gauge modes creep in. The mapping of synchronous 
perturbations onto any gauge with fixed coordinates is unambiguous as we 
have seen. However, since it is a many-to-one operation, its inverse leaves 
additional gauge freedom. Notice from equation (7), knowledge of the total 
matter or Newtonian gauge perturbations only fixes a certain combination 
of density and velocity perturbations in synchronous gauge. To remove the 
ambiguity, one must fix the initial synchronous velocity. Since an initial ve- 
locity decays with the expansion as a~^, equation (7) tells us that a different 
choice will alter the behavior of densities by an additional term proportional 
to d/a^ or during radiation domination and a~^/^ during matter dom- 
ination. This is a gauge mode. 

In summary, the "gauge invariant" approach does nothing to solve the 
gauge ambiguity; no solution is necessary since gauge choice poses no funda- 
mental problems. Yet even though there is nothing particularly deep about 
the "gauge invariant" program, it is often useful. It allows us to borrow 
Newtonian and rest frame concepts for use in any gauge. In addition, by 
providing a systematic way of mapping perturbations in a specific gauge 
onto an arbitrary gauge, one automatically determines whether the coor- 
dinates have been completely fixed, i.e. whether gauge modes have been 
entirely eliminated. 



Concepts in CMB Anisotropy Formation 9 
2.3 Newtonian Gauge Equations 

Now we are ready to state the evolution equations in an explicit form. 
To avoid unnecessary confusion, we will stay in the pure Newtonian gauge 
throughout this treatment. We will hereafter drop the superscript N with 
the understanding that all perturbation variables arc in the Newtonian 
gauge unless otherwise specified. The continuity and Euler equations for 
the photon temperature in flat space are 

6*0 = -^01-^, 

0, = k[0o + ^n^] - f (01 - Vb), 

where the term proportional to the Compton differential optical depth r 
comes from momentum conservation in the scattering. Recall that Oq = 
is the isotropic temperature fluctuation, and Oi — is the amplitude 
of the photon dipole or bulk velocity. One can see that when the optical 
depth to scattering is high, the photons become isotropic in the electron- 
baryon rest frame, i.e. the dipole moment Oi = Vb. The anisotropic stress 
of the photons 77^ is directly proportional to its quadrupole moment. Since 
scattering makes the photons isotropic in the baryon rest frame, the photon 
anisotropic stress is negligible before recombination. These equations can 
also be directly obtained from kinetic theory. In fact, they are contained in 
the Legendre decomposition of the Boltzmann equation (2) in Fourier space 
which also gives the evolution of through a hierarchy of higher angular 
moments (Wilson & Silk 1981). 

Aside from the usual velocity divergence source in the continuity equa- 
tion, there is a term dependent on the metric. This is due to the gravi- 
tational redshift effects of time dilation. As the form of the metric gij = 
—a'^dij{l + 2$e^^'^) implies, it is entirely analogous to the cosmological red- 
shift. Heuristically, the presence of matter curves or stretches space taking 
the wavelength of the photon with it. In the Euler equation, the Newto- 
nian potential W acts as a source of the dipole. Gradients in the potential 
also induce gravitational blue and red shifts as the photons fall into and 
climb out of potential wells < 0). This is countered by photon pressure 
from Oq. As the temperature rises so does the pressure which opposes the 
fall of a photon into the potential well. As the photons free stream, power 
in the dipole is converted into the higher multipole moments through the 
quadrupole via the term. 

The baryon continuity and Euler equations take on a similar form 

h = -kVb - SS, 

a (9) 
Vb = —Vb + k^ + f{Oi - Vb)/R. 
a 



(8) 
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Momentum conservation in Compton scattering gives the form of the cou- 
phng. From equation (5), the effective momentum density of a general fluid 



(4p-y/3)50i = pfe^Vfe. Thus the Compton coupling for the baryons takes on 
a similar form to the protons but is of opposite sign and is altered by a 
factor of = Spb/4:p-y. 

Again, the continuity equation is modified by gravitational effects from 
the stretching of space associated with Since the density decreases as 
the length scale cubed, the total differential effect becomes 3^. In the Euler 
equation, the velocity is damped by the expansion and enhanced by infall 
into potential wells. Particle momenta scale as due to the expansion. 
For non-relativistic particles, this causes the peculiar velocity to scale sim- 
ilarly in the absence of sources. In the fuUy-relativistic case, it causes the 
temperature to decrease as a~^. 

Decoupled components such as the massless neutrinos and cold dark 
matter follow identical evolution equations save for the absence of Compton 
coupling. The metric perturbations on the other hand feel the influence of 
the total matter perturbations, prSr = '^iPi^i^ (pt +Pt)^t = ^i{Pi + 
Pi)Vi, and ptHt — ^jP^ili, where i runs through all the particle species, 
through the generalized Poisson equations. 



which arise from the time-time + time-space and traceless space-space com- 
ponents of the Einstein equations respectively. The presence of in the first 
equation represents the conversion from physical to comoving coordinates. 
Notice that if we rewrite the equation in terms of the density fluctuation 
on the total matter rest frame, the first equation would take on the famil- 
iar non-relativistic form of the Poisson equation and simplify perturbation 
calculations. The additional term represents a relativistic effect that is im- 
portant outside the horizon. When anisotropic stress ptIIt may be ignored, 
the second of equations (10) reduces to ^ = —<P as one would expect for 
the Newtonian potential. 

2.4 Gauge Tricks: Sachs- Wolfe Example 

Examining the photon conservation equations (8) and the Poisson equations 
(10), we see that the former are simpler in the Newtonian gauge, whereas 
the latter are simpler in the total matter gauge. Let us see how gauge tricks 
developed 

in §2.2 can help us understand their joint evolution. In the matter dom- 
inated epoch, the Poisson and continuity equations reduce to 




conservation implies 



k^^ = 47TGa^pT[ST + 3-(l +PT/pT)VT/k], 

a 

k'^[^ + ^) = -STTGa^prnT, 



(10) 
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^ = 



{krj) 



5rp = 



-kVr 



as one expects from a non-relativistic analysis. Since S'^ is the density per- 
turbation on the matter rest frame, its evolution in the growing mode also 
satisfies the non-relativistic relation (xrf . Hence Vt — — 25y//c?7. Recall- 
ing its relation to the Newtonian density perturbation, we obtain 



5t = 5n 



KT] 



1 + 



12 



12 



{kriy 



2^, 



{kr]y 
/c?7 < 1 



5'^ 



For 
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adiabatic fluctuations, the intrinsic temperatiuc fluctuation, Oq = 
\5t = 1^- We shall sec that after a gravitational redshift of \P, 



this implies that the effective temperature Oo + ^ = + ^ — which is 
the famous Sachs- Wolfe (1968) result. 



3 Primary Anisotropies and the Tight CoupHng 
Approximation 

Before recombination, the Compton scattering time was so short that the 
photons and baryons behaved as a single fluid. This allows us to greatly 
simplify the treatment of anisotropy formation. Specifically, since the mean 
free path is much shorter than a wavelength of the fluctuation, the optical 
depth through a wavelength ~ f/k is large. Thus the evolution equations 
may be expanded in k/f. Employing the baryon Euler equation (9), we may 
eliminate the baryon velocity from the photon evolution equation to obtain 
the first order equation (Peebles & Yu 1970; Hu Sz Sugiyama 1995a) 

1^(1 + ii)e„4e„ = 4(1 + + (11) 

where we have dropped the higher order correction 77^ = 0{k/f) (see §3.2). 
Conceptually, this equation reads: the change in momentum of the photon- 
baryon fiuid is determined by a competition between the pressure restoring 
and the gravitational driving forces. If we ignore the time dependence of 
the baryon-photon momentum ratio R from the expansion, this equation 
describes a forced harmonic oscillator. Since scattering requires the bulk 
velocities of the photons and baryons to be equal, the effective dimensionless 
mass of the fiuid is given by mgff = 1 -|- i? to account for the inertia of 
the baryons. Baryons also contribute gravitational mass to the system as is 
evident in the appearance of mefr in the infall and dilation terms on the right 
hand side. They do not however contribute significantly to the pressure or 
restoring force of the system. 
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3.1 Acoustic Oscillations 

As an instructive first approximation, let us ignore time variations in the 
potentials and ^ and also the baryon-photon momentum ratio R com- 
pared with changes at the oscillation frequency to = kcs, where the sound 
speed (Doroshkevich, Zel'dovich & Sunyaev 1978) 

_ 1 

Equation (11) then reduces to the familiar form 

This is a simple harmonic oscillator under the constant acceleration provided 
by gravitational infall and can immediately be solved as 

Ooiv) = [0o(O) + (1 + R)^] cos{krs) + r^Oo{0) sin(A;r,) - (1 + R)W, (12) 

KCg 

where the sound horizon rg = J Csdrj ^ CgT]- The two initial conditions 
6*0(0) and ©0(0) govern the form of the acoustic oscillation. We shall see 
below that they represent the adiabatic and isocurvature modes respectively. 
Equation (12) also implies through the photon continuity equation (8) that 

01 (?y) = 3[0o(O) + (1 + R)^]cs sm{krs) + 3fc-^6>o(0) cos(/cr,). (13) 

In equations (12) and (13), lie the main acoustic and redshift effects which 
dominate primary anisotropy formation. 

3.1.1 Gravitational Infall and Redshift 

In the early universe, photons dominate the fluid and — > 0. In this limit, 
the oscillation takes on an even simpler form. For the adiabatic mode, 
6*0(0) = and 6'o(r7) = [6*0(0) + W]cos{krs) - ^. This represents an os- 
cillator with a zero point which has been displaced by gravity. The zero 
point represents the state at which gravity and pressure are balanced. The 
displacement > yields hotter photons in the potential well since grav- 
itational infall not only increases the number density of photons but also 
their energy through gravitational blueshifts. 

However, photons also suffer a gravitational redshift from climbing out 
of the potential well after last scattering. This precisely cancels the 
blueshift. Thus, the effective temperature perturbation is Oo{r]) + ^ — 
[Oo{0) + ^] cos(/crs) (see Fig. 2). The phase of the oscillation at last scat- 
tering determines the effective fluctuation. Since the oscillation frequency 
u> — kcs, the critical wavenumber k = 7T/rs{rj*) ~ tt/csT]^ is essentially at 
the scale of the sound horizon. If there is a spectrum of /c-modes, there will 
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(a) Photons: Cs=l/V3 




Fig. 2. Acoustic oscillations. Fluid compression through gravitational infall is 
resisted by photon pressure setting up acoustic oscillations. Without baryons the 
photon blueshift on infall is equal and opposite to the redshift upon leaving the 
potential leading to symmetric oscillations around zero with a Doppler effect of 
equal magnitude but 90 degrees out of phase. Baryons increase the mass of the 
fluid causing more infall and a net zero point displacement. Temperature crests 
(compression) are enhanced over troughs (rarefaction) and velocity contributions. 



be a harmonic series of temperature fluctuation peaks with km = ttitx /rsirj^) 
for the mth peak. Odd peaks thus represent the compression phase (tem- 
perature crests), whereas even peaks represent the rarefaction phase (tem- 
perature troughs), inside the potential weUs. 

As first calculated by Sachs & Wolfe (1967), the effective temperature in 
the matter dominated limit goes to ©o + = In the Newtonian frame, 
this is a direct consequence of the adiabatic initial conditions. In general, if 
kr] <^ 1 equation (8) implies 



[00 + ^] iv*) = [Oo + ^] (0) + [^~$] 



In a full calculation, the small variation in the potential at equality due to 
the change in the equation of state brings the effective temperature from 
to As this does not affect the qualitative picture, we will encorporate 
this effect as [Oq + ^]{0) = and - <?) = 0. Notice that if 6*0(0) = 
iZ/(0) = -^(0) = 0, then [0o + ^](^*) = ^^(??*) -<?(??*) ~ 2^{r]^) which yields 
a factor of 6 enhancement for isocurvature models. An alternate way of 
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Fig. 3. Baryon Drag. Baryons drag the photons into potential wells leading to a 
zero point displacement of and alternating peak heights in the rms fluctuation 
which provide a measure of Q\,h? . The fall off near o* of the amplitude is due to 
diffusion damping. 



deriving these results is to employ gauge tricks as shown in §2.4. The Sachs- 
Wolfe effect is a combination of an intrinsic temperature and a gravitational 
redshift. Since the photon density and thus the intrinsic temperature is a 
gauge dependent concept, this breakdown, but not the observable result, 
will also depend on gauge. 



3.1.2 Baryon Drag 

Though effectively pressureless, the baryons still contribute to the inertial 
and gravitational mass of the fluid meff = 1 + i?. This decreases the sound 
speed and changes the balance of pressure and gravity. Gravitational infall 
now leads to greater compression of the fluid in a potential well, i. e. a fur- 
ther displacement of the oscillation zero point (see Fig. 2). Since the redshift 
is not affected by the baryon content, this relative shift remains after last 
scattering to enhance all peaks from compression over those from rarefac- 
tion. If the baryon- photon ratio R were constant, the effective temperature 
perturbation would become 

0{n) + <f = V(l + m) cos{krs) - R^, 

o 

with compressional peaks a factor of (1-1-6-R) over the -R = case. In reality, 
the effect is reduced since i? ^ at early times (see Fig. 3). 

The evolution of the effective mass has another effect on its own. In 
classical mechanics, the ratio of energy ^meffU^A'^ to frequency lj of an 
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; Analytic Decomposition 
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Fig. 4. Analytic decomposition (i7o = 1 scale invariant adiabatic model). The ef- 
fective temperature after gravitational redshift dominates the primary anisotropy. 

Peak heights are enhanced and modulated by baryon-photon ratio R as well as 
experience a boost crossing the equality scale at £ « \/2>{2QqHq / a^qY^^rjo « 400. 
The Doppler effect is smaller and out of phase with the temperature. The ISW 
effect due to potential decay after recombination is small here but can be signif- 
icant for low matter content universes. Diffusion damping cuts off the acoustic 
spectrum at small scales. 

oscillator is an adiabatic invariant. Thus for the slow changes in mefi oc c^^, 
the amplitude of the oscillation varies as A oc oc (1 -|- R)~^/^ since 
uj (X Cs- The fundamental solutions of the oscillator equation are modified 
to be (1 + i?)-V^cos(A;rs) and (1 + i^)-^/^ sin(/crs). 

3.1.3 Doppler Effect 

Since the turning points are at the extrema, the fluid velocity oscillates 
90 degrees out of phase with the density (see Fig. 2). Its motion relative 
to the observer causes a Doppler shift. Whereas the observer velocity cre- 
ates a pure dipole anisotropy on the sky, the fluid velocity causes a spatial 
temperature variation Oi/\/2> on the last scattering surface from its line of 
sight component. For a photon-dominated Cg ~ l/-\/3 fluid, equation (13) 
tells us the velocity contribution is equal in amplitude to the density ef- 
fect. This photon-intrinsic Doppler shift should be distinguished from the 
electron-induced Doppler shift of reionized scenarios. 

The addition of baryons significantly changes the relative velocity con- 
tribution. As the effective mass increases, conservation of energy requires 
that the velocity decreases for the same initial temperature displacement. 
Thus the relative amplitude of the velocity scales as Cg ■ In the toy model of 
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a constant baryon-photon density ratio R, the oscillation becomes Oi/y/S = 
^^^(l + 3i?)(l + sin(/crs). Notice that velocity oscillations are sym- 

metric around zero unlike the temperature ones. Thus compressional peaks 
will rise clearly above the velocity oscillations if R is large. Even in a uni- 
verse with Qbh'^ given by nucleosynthesis, R is sufficiently large to make 
velocity contributions subdominant (see Fig. 4). 

3.1.4 Acoustic Driving Effects and Isocurvature Models 

Whenever the non-relativistic matter is not the dominant dynamical com- 
ponent, the potentials ^ and ^ become time-dependent. For example, when 
the universe is radiation dominated, pressure and entropy alter the behav- 
ior of the gravitational potential. External sources from topological defects 
can also play a role. The effects of potential evolution can be separated into 
those that occur before last scattering and thus affect the acoustic oscilla- 
tions and those that occur afterwards which affect the gravitational redshift 
of the photons along their free-streaming trajectories. We will defer consid- 
eration of the latter to §4 where we consider effects between recombination 
and the present. 

As a simple example of the driving effects of potential evolution, let us 
again consider adiabatic fluctuations where the initial temperature Oo{0) 
and curvature fluctuations ^(0) are finite constants related by the Poisson 
equation and 00(0) = 0. Inside the sound horizon, pressure prevents grav- 
itational infall in the photon-baryon and neutrino systems. Unless CDM 
dominates, energy density perturbations 6pT decay with the expansion and 
consequently can no longer maintain a constant gravitational potential. 
Counterintuitively, this decaying potential can actually enhance tempera- 
ture fluctuations through its near resonant driving force. Since the potential 
decays after sound horizon crossing, it drives the first compression without 
a counterbalancing effect on the subsequent rarefaction stage (see Fig. 5). 

Furthermore, the dilation effect from a decaying space curvature ^ also 
enhances the acoustic fiuctuations. Heuristically, the overdensities which es- 
tablish the potential well "stretch" the space-time fabric. As the potential 
well decays, it re-contracts. Photons which arc caught in this contraction 
find their wavelength similarly contracted, i.e. blueshifted. Thus a differen- 
tial change in $ leads to a dilation effect, Oq — Combined, the driving 
effects of infall and dilation yield an effect of order — ^ + <P ^ 
or 5 times the Sachs- Wolfe effect when the perturbation crosses the sound 
horizon (see Fig. 4,5 and Hu & Sugiyama 1995c). Since this effect only oc- 
curs for modes which cross before CDM domination, the amplitude of this 
boost is sensitive to the matter-radiation ratio of the universe. In particular, 
lowering f2oh'^, raising the massless neutrino number or energy density, or 
making some of the cold dark matter relativistic at early times through eV 
mass neutrinos (see Fig. 6b; Ma & Bertschinger 1995; Dodelson, Gates & 
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Fig. 5. Driving effects. The time evolution of the potential can enhance the am- 
plitude of the acoustic oscillation by driving effects. In the adiabatic case [(a) solid 
lines, numerical results and (b) heuristic picture], the potential enhances the first 
compression through infall and decays leaving the oscillator strongly displaced 
from the zero point. In the (baryon) isocurvature case [(a) dashed line], the po- 
tential grows from zero and stimulates a sine mode. The first extrema here is 
suppressed due to the fact that the gravitation driving begins near sound horizon 
crossing. Notice that for the first cycle of the adiabatic and isocurvature oscilla- 
tions, the gravitational force mimics a driving force of approximately twice the 
natural period. 

Stebbins 1995), can move the boost to larger scales and thus enhance the 
first few acoustic peaks. 

Similar resonant effects can occur in other situations. For example, many 
alternate scenarios follow a typical isocurvature behavior. The hallmark of 
isocurvature type models is that curvature fluctuations $ are zero or at least 
suppressed outside the horizon. The curvature fluctuation grows by causal 
processes as the fluctuation crosses the horizon. This isocurvature behavior 
occurs for entropy fluctuations where matter and radiation perturbations 
are counterbalanced in the initial state, whether through the baryons (Hu 
& Sugiyama 1995b), hot dark matter (de Laix & Scherrer 1995) or axions 
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Fig. 6. Heights of the peaks in a scale invariant adiabatic model. The heights of 
the peaks is determined by the baryon drag and acoustic driving effects and so 
are sensitive to Q^h^ and Qoh? . Baryons increase compression in potential wells 
causing the peak heights to alternate. Lower matter content causes more potential 
decay driving the oscillation. 



(Sugiyama & Gouda 1992), texture models (Crittenden & Turok 1995, Dur- 
rer, Gangui & Sakellariadou 1995), and is at least part of the story for string 
models (Hindmarsh, private communication, but see also Albrecht, Coulson 
& Ferreira 1995). The qualitative effect on isocurvature conditions is easy to 
see. Since the potential grows from zero to a maximum near sound horizon 
crossing and then decays due to radiation pressure, the force drives the sine 
harmonic of oscillations (see Fig. 5a). The result is that isocurvature acous- 
tic peaks are 90 degrees out of phase with their adiabatic counterparts. The 
first peak occurs at larger scales that the corresponding adiabatic peak but 
generally has a suppressed amplitude since the corresponding mode is just 
at the sound horizon krg = 7r/2 and has not experienced the full effect of the 
driving term. These features may serve to distinguish the two basic types 
of scenarios in a manner independent of the details of the given model. 

Notice that in these examples, oscillations give a coherent spectrum of 
peaks in integral or half integral multiples of J udrj = krs{r]^). Since each 
k mode evolves independently, this coherence in the temporal phase is the 
result of a special timing in the gravitational impulse. In all the examples 
considered above, the gravitational force kicks in between horizon cross- 
ing and sound horizon crossing. This explains the coherence in k and the 
fundamental scale in the problem (?]*). More exotic scenarios may not 
preserve this coherence. If the source of the gravitational potential has a 
more random temporal behavior, the sine and cosine acoustic modes will 
be stimulated incoherently with equal likelihood. Combined, the two modes 
would lead to a smoothed out rms temperature fluctuation below the sound 
horizon. Preliminary calculations indicate that this may be case for cosmic 
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string models due to complicated behavior of the string network inside the 
horizon (Albrecht et al. 1995; Magueijo et al. 1995). 

These considerations can be made quantitative by solving the oscilla- 
tor equation under the influence of an arbitrary but known gravitational 
forcing function either analytically (Hu & Sugiyama 1995a) or numerically 
by modeling the gravitational source (e.g. Seljak 1994; Crittenden & Turok 
1995). Simple closed form solutions for standard adiabatic and isocurvature 
models are presented in Hu & Sugiyama (1995c). 

3.2 Diffusion Damping 

In reality, the photons and baryons are not perfectly coupled since the pho- 
tons possess a mean free path in the baryons Ac ~ due to Compton 
scattering. As the photons random walk through the baryons, hot spots and 
cold spots are mixed. Fluctuations thereafter remain only in the unscattered 
fraction causing a near exponential decrease in amplitude as the diffusion 
length ~ y/NXc = Vv^c or kri ^ \/'^/v overtakes the wavelength. 

To be more specific, diffusion causes heat conduction and generates vis- 
cosity in the fiuid (Weinberg 1972). As photons from regions of different 
temperature meet, anisotropics form leading to a quadrupole moment or 
anisotropic stress in the fluid. If the diffusion length is well under the hori- 
zon, it overtakes the wavelength of the fluctuation when f/k = krj ^ 1. 
Thus the optical depth through a wavelength is still high and the pertur- 
bative expansion of §3.1 still holds. This fact allows us to extend the tight 
coupling approximation for acoustic modes to last scattering. For small scale 
modes in which the optical depth through a wavelength f/k is not high at 
recombination, all acoustic oscillations will already have damped away (Hu 
& Sugiyama 1995c). 

Photon diffusion is a second order effect and damps acoustic oscillations 
as exp[— (A;/A;d)^] with the damping wavenumber 



where /2 accounts for the subtle effects in the generation of anisotropic 
stress (see Fig. 4). Kaiser (1983) showed that /2 = 9/10 due to the angular 
dependence and 3/4 if the additional effects of polarization are included. 
The two processes aid the generation of the quadrupole moment and hence 
increase the viscous damping of the acoustic oscillations. 

The baryons are dynamically coupled to the photons by momentum 
exchange in Compton scattering from fd = f/Rin equation (9). If fd//c 3> 1, 
the baryons will be dragged in and out of potential wells with the photons. 
This process destroys the baryonic acoustic oscillations as well and is known 
as Silk (1968) damping. Since coupling requires coevolution in the number 
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density 5b = = SOq, only entropy fluctuations S = 5}, — ^5^ survive 
diffusion damping. 

Notice also that the ionization history factors in the diffusion length 
through the mean free path f~^. As we shall now see, at recombination the 
mean free path and hence the diffusion length increases substantially but 
does not approach the horizon scale. 

3.3 Decoupling 

The CMB anisotropy today is simply the acoustic fluctuation at last scat- 
tering, modified by diffusion damping and free streamed to the present: 

[0 + ^]{r]o,k,ii) ^ [Oo + ^ -il^Oi]{T]„k)V^{k)e'''i'^''*-'^°\ (14) 

for flat space, where k • 7 = k/i, 

(■no 

V^{k)= dryV^e-l'^/'^^^^)] , (15) 
Jo 

is the diffusion damping factor and the visibility function = fe~'^ is the 
probability of last scattering within dr] of 77. Its peak is near t{z^) = 1. The 
damping factor is dependent only on the background cosmology and fits 
across the whole range of parameter space are presented in Hu & Sugiyama 
(1995c). We shall examine the meaning and implications of each of these 
terms below. 

The baryons do not decouple from the photons precisely at last scatter- 
ing. As we have seen in §3.2, the coupling strength is altered by a factor of 
R, fd = f / R. By analogy to the photon case, we can define a drag optical 
depth Td and the end of the drag epoch as Td{zd) = 1 (Hu & Sugiyama 
1995c). After this point, Compton drag on the baryons can no longer pre- 
vent the gravitational infall of the baryons. Because i?(2;*) ~ 0.3 for standard 
recombination and a big bang nucleosynthesis value for Q^h?, and Zd ap- 
proximately coincide for the standard case. However for reionized scenarios 
R{z^) ^ 1 and the baryons decouple from the photons long before last scat- 
tering. It is thus no longer appropriate to consider the photons and baryons 
as tightly coupled at last scattering. We will develop new techniques to han- 
dle this situation in §4. From the drag depth r^, a drag visibility function 
Vb ~ fde~'^'^ can be constructed (see Hu & Sugiyama 1995c for the small 
correction due to expansion damping). The final scale for Silk damping is 
obtained from the visibility function and the diffusion length as in equation 
(15) for the photons. 
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Fig. 7. Location of the peaks. The location of the peaks is determined by the 
angle the sound horizon subtends at last scattering. In an open universe, geodesic 
deviation makes this angle much smaller than in a flat universe. As this angle 
involves a ratio of the sound horizon to the (angular diameter) distance to last 
scattering, distance scale changes through A and h have little effect on its angle. 
The drop in the lowest multipoles for the open Qq = 0.1 is due to the lack 
of supercurvature scale power in the otherwise scale invariant initial conditions 
assumed here. Models here have h = 0.5 and Qb = 0.05. 

3.4 Projection Effects and the Anisotropy Spectrum 

The presence of the exp{ikiJ,A'r]) term in equation (14) represents the free 
streaming of the photons in flat space. Photons travel at the speed of light so 
that the number of wavelengths traveled between last scattering and today is 
kAi]. However, the phase of the wavefront changes only in the perpendicular 
direction so that the observed phase change will depend on the line of sight 
(f) = kArjii. Since inhomogeneities at the last scattering surface appear as 
anisotropies on the sky today, this can equivalently be viewed as a simple 
projection of the plane wave on the sphere. The anisotropy is expressed 
through the decomposition of the plane wave into multipole moments i ~ 
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Fig. 8. Parameter sensitivity. A scliematic representation based on a scale invari- 
ant adiabatic scalar model. Features in an open model are shifted to significantly 
smaller angles compared with A and i7o = 1 models, represented here as a shift 
in the i axis. Isocurvature models behave similarly with respect to the acoustic 
oscillations save that the peaks are 90 degrees out of phase. The spectrum may 
also be tilted by changing the initial power spectrum from the scale invariant 
assumption. 



-1 



= [[Oo + m*.k) + 01 (ry,, k)^^'^ V^ik)MkArj). (16) 

Since ji peaks at £ ~ kArj, it is clear that free streaming just projects 
the physical scale /c ~ onto an angular scale as 6* ~ X/Arj. The total 
power in the ith multipole is obtained by integration over k modes, Cg — 
(2/71) Jk^\ei{r]Q, A;)|V(2£ + l)^dlnk . 

The physical content of the generalization to open universes is now ob- 
vious. Since free streaming merely represents a projection, one replaces 
the comoving distance Ar] with the comoving angular diameter distance 
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re = \K\-^/'^ siBh[\K\^/^Ar]], where the curvature K = -H^{1 - - r2yi). 
For a, K > closed geometry, merely replace the sinh with sin. Notice that 
in the closed case, the location of the first peak has oscillates with the curva- 
ture (White & Scott f995). Putting this together, the acoustic peaks occur 
at 



Formally, this arises by replacing jj^ with the radial eigenfunction of an 
open universe (Wilson 1983). Notice that the peak locations only depend 
on the background cosmology once the adiabatic or isocurvature nature of 
the fluctuations has been established. Indeed, the ratio of the peak locations 
themselves can be employed to separate the two (Hu & Sugiyama 1995b). 
The dominant factor in the peak locations is the curvature of the universe 
which can make the same physical scale subtend a much smaller angle on the 
sky (see Fig. 7). Here we have a classical cosmological test of the curvature: 
knowing the physical scale (sound horizon at last scattering) and the redshift 
{z^ f» 1000) of the acoustic peak, we measure its angular extent. The inferred 
angular diameter distance tells us the curvature. If the peak location is 
known to sufficient accuracy, we can also measure Qi^h?, f2oh^ and A through 
their effects on the sound horizon and angular diameter distance (see Figs. 
6 & 7). Analogous information can be obtained from the location of the 
diffusion damping scale even if oscillations are not apparent due to random 
driving effects. 

The heights of the peaks also contain important, if less model inde- 
pendent information. As we have seen, every other peak is boosted by the 
baryon content due to baryon drag and scales inside the horizon at equal- 
ity probe f^o^^ or more generally the matter-radiation ratio which includes 
information on the neutrino content and mass. Since these two effects have 
relatively robust features, alternating peak heights and a strong height boost 
at matter-radiation equality, there is hope of obtaining information out of 
the heights of the peaks as well as their location for a wide class of models. 
As an example, we present a schematic picture for the anisotropy spectrum 
in Fig. 8 based on the popular scale invariant adiabatic model. Note that 
changing the overall dynamics from Qq = 1 through flat A models to open 
models is similar to shifting the spectrum in angular space toward smaller 
angles. 




jTT, adiabatic 
{j — 1/2)tt. isocurvature 
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4 Secondary Anisotropies and the Weak Coupling 
Approximation 



We have been assuming up to this point that CMB fluctuations are frozen 
in at a redshift ~ 1000. Several processes in the foreground of recombi- 
nation could alter the anisotropy spectrum. These are generally known as 
secondary anisotropies. Just as the tight coupling approximation assisted in 
the calculation and interpretation of primary anisotropies, the weak cou- 
pling approximation helps in understanding secondary anisotropies (Hu & 
White 1995). 

In abstract form, anisotropy generation is governed by internal sources 
from photon fluctuations at last scattering Sprim ^ind sources external to 
the photon system in the foreground of recombination Ssec by 



in flat space with an appropriate generalization of the radial eigenfunction 
for an open geometry (Wilson 1983). Equation (17) just states that we ob- 
serve the projection of the source at ry on a shell at a distance At] — rjo —rj. 
Whether the sources or the phase of the wave varies more rapidly distin- 
guishes the tight from the weak coupling regime. For primary anisotropies, 
the source is localized over a short range of time around recombination. 
Since f//c ^ 1, ji can be taken out of the integral (17) leading to the simple 
projection of equation (16). 

Secondary anisotropies do not necessarily possess this property. Most 
cosmological effects aside from decoupling take on the order of an expansion 
time at the relevant epoch to be completed. In these cases, Ssec may be taken 
to be slowly varying at small scales and removed from the integral in (17). 
Since 



£ more rapidly than a simple projection of the source would imply. This 
just reflects the fact that it contributes across many wavelengths of the 
fluctuation allowing contributions from crests and troughs to cancel. We 
will now discuss the main sources of secondary anisotropies Ssec- 




(17) 




(18) 



in the weak coupling limit where 
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4.1 Gravitational Effects 



Even in the absence of reionization, secondary anisotropics can be gener- 
ated by gravitational redshift effects between recombination and today. The 
differential redshift from ^ and dilation from ^ discussed above must be 
integrated along the trajectory of the photons Sgec = ^ — Tensor fluctua- 
tions can also give rise to anisotropics through this mechanism. We thus call 
the combination the integrated Sachs-Wolfe (ISW) effect. Notice that these 
effects only occur if the metric fluctuation is time varying. We can separate 
this general mechanism for anisotropy formation into categories based on 
the reason for the time evolution. There are four possibilities to consider: 
time evolution due to the radiation content (early ISW effect), due to the 
expansion (late ISW effect), due to non- linear evolution (Rees-Sciama ef- 
fect) and any more exotic sources such as gravity waves or defects (sourced 
ISW). Additionally, beyond linear theory gravitational lensing may affect 
the CMB. Since this merely shuffles power in anisotropics around scales at 
the arcminute level for a CDM type model, we will not further consider 
the effect and refer the interested reader to Seljak (1995a) and references 
therein. 



4.1.1 Early ISW Effect 



The early ISW effect is the direct analogue of the acoustic driving effect 
of §3.1.4 except that the photons are in the free streaming rather than 
the tight coupling regime. For adiabatic conditions, the potential decays 
after horizon crossing in the radiation dominated limit. For isocurvature 
conditions, it grows outside the horizon and then decays as in the adiabatic 
case. In general, this effect will smoothly match onto the acoustic peaks 
due to similarities in the cause of their generation. Due to the later time of 
generation, this effect influences larger scales than the acoustic peak but is 
cut off above the equality scale. Furthermore, since it arises from a distance 
closer to ourselves, the same physical scale subtends a larger angle on the 
sky. Together these considerations imply that the early ISW effect flUs in 
the anisotropy on scales just larger than the flrst acoustic peak (see Fig. 8) . 
It serves to broaden the rise and shifts the location of the flrst peak in 
the spectrum to larger scales. Unfortunately, this effect satisfies neither the 
tight nor the weak coupling approximation since the decay time ~ rj and 
the wavelength ~ are by definition comparable for this horizon crossing 
effect. 
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4.1.2 Late ISW Effect 

If the universe has a non- vanishing curvature or cosmological constant, even- 
tually these will dominate the expansion rate. Under the rapid expansion, 
density fluctuation 5pT decays taking the potential with it. The decay takes 
on the order of an expansion time at the end of matter domination indepen- 
dent of the wavelength. Since the photons free stream, they travel across 
many wavelengths of the perturbation on scales smaller than the horizon. 
If the potential decays while the photon is in an underdense region, it will 
suffer an effective redshift rather than a blueshift. Contributions from over- 
dense and underdense regions will cancel and damp the ISW effect on small 
scales. This is the hallmark of the weak coupling regime and is mathemati- 
cally expressed through cancellation in the integral (18). 

For a fixed Qq, the decay epoch occurs much later in flat Qa + ^^o = 1 
models than open ones. Thus A models will suffer cancellation of late ISW 
contributions at a much larger scale than open models. In fact, for reasonable 
Qq>Q.1 decay the decay has already started at the quadrupole (see Fig. 7,8; 
Kofman & Starobinskii 1985; Hu & White 1995). In summary, the epoch 
that the universe exits the matter dominated phase is imprinted on the 
CMB by the late ISW effect. 

4.1.3 Rees-Sciama Effect 

Even for sjv = 1 CDM dominated universe, potentials only remain con- 
stant in linear perturbation theory. The second order contribution has been 
shown to be negligibly small (Martinez-Gonzalez, Sanz & Silk 1992). The 
effect in the non-linear regime has been estimated with N-body simulations 
through power spectrum techniques (Seljak 1995) and ray tracing techniques 
(Tuluie, Laguna & Anninos 1995). The general conclusion is that this effect 
does not become comparable to the primary signal until well into the diffu- 
sion damping tail. In the absence of reionization, it should therefore present 
no problem for the extraction of main features such as the curvature and 
the baryon content from the acoustic peaks. It may however complicate the 
extraction of cosmological parameters from more subtle effects, for example 
the neutrino mass. 

4.1.4 Sourced ISW Effect 

This catch-all category contains all other gravitational redshift effects on 
the CMB. Tensor perturbations in the metric, i.e. gravity waves, give rise 
to dilation effects as the photons free stream. Due to the nature of the 
metric distortion, they leave a quadrupole signature in the CMB which 
thereafter is projected onto higher multipoles. Combined with cancellation 
effects, this implies that the tensor spectrum thus typically exhibits a sharp 
drop in power from the quadrupole. Since contributions only arise between 
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Fig. 9. Gravitational Wave ISW Effect. Tlie metric distortion caused by a grav- 
itational wave induces a redshift/dilation effect on CMB photons. The spectrum 
is cut off below the angle subtended by the horizon at last scattering since gravity 
waves only affect free streaming photons significantly. Due to the signature of the 
metric distortion, it originates as a source of the quadrupole anisotropy. 



last scattering and the present, there is also a cut off at the angle the horizon 
subtends at last scattering. Below this scale, the gravity waves have already 
redshifted away by last scattering. Detailed calculations of the spectrum 
were first carried out by Crittenden et al. (1994) and the spectrum shown 
in Fig. 9 is from Hu et al. (1995). For inflationary models, there exists a 
consistency relationship between the amplitudes of the scalar and tensor 
modes and the slope of the power spectrum (see e.g. Steinhardt 1995). 
There is hope that detailed measurements of the CMB spectrum can thus 
provide a strong test of the inflationary scenario (see Lindsey et al. 1995 
and references therein). 

Cosmological defects may also act as an external source of gravitational 
potential perturbations. Pen, Spergel & Turok (1993) find that in such mod- 
els a scale invariant form for the anisotropy results from this effect and can 
mimic the infiationary prediction. As in the case of the late ISW effect, the 
magnitude drops from 2^ to zero as the characteristic time scale for the 
change in ^ becomes longer than the light travel time across a wavelength. 
For defect models Pen, Spergel & Turok (1993) estimate the amplitude as 
approximately Since this effect occurs around horizon crossing for a 
given mode due to the isocurvature conditions, the largest spatial modes 
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Fig. 10. Reionization damping, Doppler, and Vishniac eflFects (scale invariant adi- 
abatic i?o = 1, i^b = 0.05, h = 0.5, COBE normalized model). As the ionization 

level between recombination and the present increases, primary anisotropics are re- 
duced by diffusion and rescattering. For sufficiently high ionization, here Zi > 100, 
the Doppler effect on the new last scattering surface can regenerate some fluctua- 
tions at intermediate scales. The Vishniac effect dominates at very small angular 
scales. Due to its second order nature, it contributes even if the ionization redshift 
is low. The Zi — 5, 10 curves are indistinguishable for the primary anisotropics as 
are the Zi = 100, 1000 for the Vishniac effect. 

have not been enhanced and can lead to a drop in power at the lowest 
multipoles. 



4.2 Scattering Effects 

An early round of structure formation may be able to reionize all or part of 
the universe at high redshift. If this occurs at 

z>z^^ 102(^2o/^V0•25)^/^(a;e^26/iV0•0125)"^/^ 

where r(z*) — 1 then CMB anisotropics will be drastically influenced 
by scattering. The model-independent information stored in the primary 
anisotropy spectrum will be greatly reduced. It is replaced by more detailed 
information and clues about the evolution of structure in the universe. 
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4.2.1 Reionization Damping 

The most dramatic effect of reionization is the increase in the photon dif- 
fusion length. Since last scattering is delayed until the Compton mean free 
path approaches the horizon (or equivalently when the scattering rate f 
drops below the expansion rate d/a), the diffusion length at last scattering 
is the horizon scale ~ ~ ^Jrja/a ~ rj. Intrinsic photon fluctua- 

tions such as the acoustic oscillations will be damped by diffusion below the 
horizon scale. Moreover, since acoustic oscillations appear only below the 
sound horizon no oscillations will be apparent in the reionized spectrum. 

If reionization is not sufficiently early, some trace of the acoustic oscilla- 
tions may remain in the spectrum. Recall that diffusion damping works since 
rescattering of photons arriving from directions with different intrinsic tem- 
peratures varies destroys the anisotropy. Fluctuations are only retained in 
the unscattered fraction e~'^ . If the total optical depth between recombina- 
tion and the present is r < 1, primary and secondary scattering anisotropies 
may be present in the spectrum (see Fig. 10). In this case, the information 
contained in the CMB would truly be enormous but also difficult to extract. 



4.2.2 Cancelled Doppler Effect 



After the Compton drag epoch j^z^ lQb{QQh'^)^^^Xe ^^'^ where Td{zd) = 1 
(Hu & Sugiyama 1995c), baryonic gravitational instability can no longer 
be prevented by the photons. Collapse of the baryon density fluctuations 
implies that baryon peculiar velocities will create a Doppler effect in the 
CMB. However the weak coupling approximation tells us that these too 
will be damped as the photon traverses many crests and troughs of the 
perturbation at last scattering. In fact, cancellation for the Doppler effect is 
particularly severe. If the perturbation wavevector is perpendicular to the 
line of sight k ± 7, cancellation is avoided. However, flows are irrotational 
in linear theory so V5 || k. Since the velocity is then perpendicular to the 
line of sight, no Doppler effect arises 7 ■ v;, = v^'^ ■ k = (see Fig. 11). 
Fluctuations only survive if there are variations in the velocity or optical 
depth through last scattering. Thus the Doppler source is related to the 
derivative of these quantities as Ssec = e~'^[Vhf + Vi,f]/k at small scales 
where recall V(,(x) = —iVt,e'^^'^'k. This leads to a suppression by a factor of 
approximately {kr]^)~^ (Kaiser 1984). Cancellation under the weak coupling 
approximation yields a further suppression via equation (18). In general, 
later last scattering implies a greater krj^, a larger cancellation scale, and 
hence smaller Doppler fluctuations. Nevertheless, in Fig. 10 the slight upturn 
at intermediate scales for the high ionization case is due to this effect. 



O 
I— ' 

O 

< 
< 



o 



03 





e velocity 
overdensity 



unscattered J 
blueshifted J 



\U;/ Vishniac 
W Mechanism 



Observer ? 



e velocity 



I : 
f / 



redshifted J 




overdensity 



blueshifted y 



% I Cancellation 



1 t 

Observer V 



Concepts in CMB Anisotropy Formation 31 



4.2.3 Vishniac Effect 

Since the first order Doppler fluctuation is severely suppressed due to "co- 
incidental" geometrical reasons other effects may dominate the anisotropy 
at small scales. Variations in the optical depth due to density fluctuations 
across the last scattering surface may alow the Doppler effect to escape 
severe cancellation. The increased probability of scattering causes a prefer- 
ential generation of Doppler fluctuations in overdense regions. The ensuing 
anisotropy is not as severely damped as the flrst order contribution. This 
is because the velocity fleld can be parallel to the line of sight while the 
density variation is perpendicular to the line of sight (see Fig. 11). This 
effect can also be calculated under the weak coupling approximation. It is 
slightly more complicated than linear effects since it entails mode coupling 
between the density and velocity flelds. An explicit expression is given in Hu 
Sz White (1995). For a CDM model, we flnd that Vishniac contributions do 
indeed exceed linear contributions at sufficiently small scales. Furthermore, 
since second order effects are strongly peaked to late times, even minimally 
ionized scenarios carry a significant fraction of the maximum signal (see Fig. 
10). 



4.2.4 Cluster Sunyaev-Zeldovich Effect 

Clusters provide a non-linear analogue of the Vishniac effect. Here the hot 
cluster provides the variation in the optical depth which causes preferential 
scattering. The Doppler effect due to the peculiar velocity of the cluster 
yields an anisotropy known as the kinematic Sunyaev-Zel'dovich (SZ) effect 
(Sunyaev & Zel'dovich 1972). For an individual cluster the temperature 
fluctuation is of order TcVc, where the optical depth through a cluster is 
of order Tc ~ 0.1 — 0.01 and the peculiar velocity Vc ~ few x 10~^. This 
provides an interesting way of measuring the peculiar velocity of a cluster 
without introducing the problems associated with determining the distance 
scale (see e.g. Sunyaev & Zel'dovich 1980, Haehnelt & Tegmark 1995). The 
average effect is much smaller however and probably does never dominates 
the anisotropy spectrum. Persi et al. (1995) estimate the effect in a CDM 
model from hydrodynamic simulations and find that this non-linear Doppler 
effect is small in comparison to the second order Doppler (Vishniac) effect 
Compton scattering off hot electrons also produces spectral distortions 
which look like anisotropics to any experiment confined to low frequency 
measurements. As discussed in §2.1, Compton scattering exchanges energy 
between the electrons and photons to order Tg/mg. This upscatters photons 
from the Rayleigh- Jeans to the Wien regime and leads to a temperature dis- 
tortion of {AT/T)rj = -2y « -2rcTe/me ^ 10"^ - lO'^ for an individual 
cluster. Again, the rms fluctuations would be much smaller. Much effort 
has been expended to estimate the fluctuations caused by this thermal SZ 
effect with varying results (Barbosa et al. 1995; Makino & Suto 1993; Cole 
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& Kaiser 1988). Recently, empirical modeling of clusters has predicted that 
the anisotropy at arcminutes should be on the order of {AT/T)jij < 10~^ 
(Ceballos & Barcons 1994). Moreover the signal is in large part due to bright 
and easily identifiable clusters. If such known clusters are removed from the 
sample, the anisotropy drops to an entirely negligible level. 

4.2.5 Inhomogeneous Reionization 

Finally, let us mention another variant of the Vishniac-Sunyaev-Zel'dovich 
mechanism. The optical depth could vary due to inhomogeneities in the 
ionization fraction last scattering. Note that inhomogeneities well 

before optical depth unity have little eff'ect on the CMB due to the cut 
off from the visibility function. Details of this effect will of course depend 
on the exact model for reionization, i.e. the extent of the inhomogeneities. 
Small inhomogeneities may be expected to behave as the Vishniac effect; 
large inhomogeneities like the kinetic SZ effect. 

5 Discussion 

It should now be clear that a great wealth of information about cosmology 
and the model for structure formation lies waiting to be observed in the 
CMB anisotropy spectrum. The location of the first acoustic peak provides 
a robust classical test for curvature in the universe if the fluctuations are 
known to be either adiabatic or isocurvature. The only caveat here is that 
if the universe suffered early reionization the acoustic effect may be so sup- 
pressed that it becomes unobservable. In this case, the CMB will place a 
lower limit on the epoch of reionization. It is very likely that the curvature 
and/or thermal history of the universe will be measured in the near future. 

Barring early reionization, once we obtain precise measurements of the 
first peak and beyond, we should be able to extract much more information. 
The relative locations of the higher peaks can cleanly separate adiabatic and 
isocurvature models independently of the curvature and thermal history. A 
precise measurement of the location of the peaks can supply information on 
A and Qoh'^. The heights of the peaks yield even more information if some 
minimal assumptions are made for the theory of structure formation. In- 
deed, requiring consistency with large scale structure measurements should 
eventually fix the model quite precisely. The relative heights of the peaks 
give a robust probe of the baryon content Q^h"^ . In any given model, the 
absolute heights yield a constraint on the matter content of the universe 
Q^h'^ and perhaps the number, temperature and mass of the neutrinos. 
Combining these pieces of information we can infer the Hubble constant. 

It would seem that all the fundamental cosmological parameters are 
encoded in the CMB anisotropy spectnnn. Yet, even discoiniting the pos- 
sibility of early reionization, how likely is it that we will precisely measure 
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them with the next generation of experiments? The results of Jungman et 
al. (1995) suggest that at least the curvature can be measured to better than 
5% accuracy with a full sky map to half a degree resolution. Precisely how 
much information can be extracted will depend in the end on how severe 
foreground contamination from synchotron radiation, free-free emission, in- 
terstellar cold dust, and radio point sources will be. Many of the effects 
described here will require AT/T to 10~^ accuracy to measure definitively. 
With sufficient frequency coverage, there is hope of distinguishing the back- 
ground from the foreground signal to employ at least the clean patches of 
the sky for cosmology. 
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